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Abstract. 

The geometrical critical behaviour of the two-dimensional Q-state Potts model is 
usually studied in terms of the Fortuin-Kasteleyn (FK) clusters, or their surrounding 
loops. In this paper we study a quite different geometrical object: the spin clusters, 
defined as connected domains where the spin takes a constant value. Unlike the 
usual loops, the domain walls separating different spin clusters can cross and branch. 
Moreover, they come in two versions, 'thin' or 'thick', depending on whether they 
separate spin clusters of different or identical colours. For these reasons their critical 
behaviour is different from, and richer than, those of FK clusters. We develop a transfer 
matrix technique enabling the formulation and numerical study of spin clusters even 
when Q is not an integer. We further identify geometrically the crossing events which 
give rise to conformal correlation functions. We study the critical behaviour both in 
the bulk, and at a boundary with free, fixed, or mixed boundary conditions. This 
leads to infinite series of fundamental critical exponents, ht 1 —t 2i 2t 1 in the bulk and 
^i+2(fi-f 2 ),i+«i a t the boundary, valid for < Q < 4, that describe the insertion 
of l\ thin and li thick domain walls. We argue that these exponents imply that the 
domain walls are 'thin' and 'thick' also in the continuum limit. A special case of the 
bulk exponents is derived analytically from a massless scattering approach. 



PACS numbers: 64.60.De 05.50+q 
1. Introduction 

Many of the key developments in the study of two-dimensional critical phenomena 
originate from the study of only a very few lattice models [TJ [2] . These include dimer 
coverings, the Ising and six- vertex models, the O(n) model, and the Q-state Potts model. 
The latter two models are particularly important, as they can be formulated in terms 
of geometrical degrees of freedom that in turn describe extended fluctuating objects, 
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such as domain walls in magnets, percolation clusters, and polymers adsorbed on walls 
and interfaces [3]. These objects are in fact conformally invariant whenever n or \[Q is 
comprised in the range —2 < n < 2, and their fluctuations are characterised by critical 
exponents. 

The purpose of this paper is to revisit the geometrical formulation of the two- 
dimensional Potts model. In particular, we shall define a set of geometrical observables 
which are different from those considered in most of the existing literature. These 
observables are simply the domain walls in the formulation of the Potts model in terms 
of Q-component spins. Not only are they in many ways more natural and physically more 
relevant than the Fortuin-Kasteleyn (FK) clusters jl] usually considered, but they also 
turn out to have a richer critical behaviour and a more complete set of critical exponents. 
These features stem in part from the fact that there are actually two different types of 
spin domain walls (thin and thick), as we shall explain shortly. 

The Q-state Potts model is defined by the partition function 

z = e n ex p i K6 w) > w 

<J (ij)&E 

where K is the coupling between spins <7j = 1, 2, . . . , Q along the edges E of some lattice 
C For simplicity we shall take £ to be the square lattice in the computations below, 
whereas we use the triangular lattice in the figures. The Kronecker delta function S aij(Tj 
equals 1 if cr, = <Tj, and otherwise. 

The usual route p[] is to write the obvious identity 

exp (K8 ait(Tj ) = 1 + vS auaj , (2) 

with v = e A — 1, and to expand Z in powers of v. The result is an expression of Z 
as a sum over FK clusters with weight v per unit length and fugacity Q per connected 
component. Alternatively one can think of the FK clusters in terms of their surrounding 
hulls, which are self and mutually avoiding loops with fugacity n = y/Q [2J. Most 
features of these FK clusters and loops, in the critical regime —2 < n < 2, are by 
now under complete control, thanks to the combined powers of Conformal Field Theory 
(CFT) and Schramm-Loewner Evolution (SLE) [5J [6] . In particular, the loops behave 
like the SLE trace in the continuum limit [B], and viewing them as contour lines of 
a (deformed) Gaussian free field leads to the Coulomb Gas (CG) approach to CFT 
[5j [7]. Our understanding of the critical properties of FK clusters and loops can be 
considered almost complete, although some of their more intricate observables — such as 
the backbone [8] and shortest-path [9] dimensions — are still unknown. 

The expansion in powers of v has some pleasing features, notably that the FK 
clusters coincide with percolation clusters in the formal limit Q — > 1. But apart from 
that it is somewhat artificial in view of the original formulation ([I]) in terms of Q- 
component spins. In particular, albeit two lattice sites belonging to the same FK cluster 
will necessarily have the same spin value, the converse is not true. It would therefore 
seem more natural to perform the expansion in powers of e K and consider as basic 
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Figure 1. A configuration of the Q = 3 Potts model, and the corresponding set of 
branching domain walls. 

geometrical observables the connected domains with a constant value of the Potts spin, 
henceforth referred to as spin clusters. 

The properties of these spin clusters have as a rule remained ill understood. This 
is particularly frustrating, since those are the very clusters that one would observe in 
an actual experiment on a magnetic alloy in the Potts universality class. The case 
of the Ising model Q = 2 is an exception to this rule [TUl ITT] , but this is due to its 
"coincidental" equivalence to the O(n) vector model with n — 1. Indeed, defining the 
Ising spins on the triangular lattice, the corresponding 0(n) model is described by self 
and mutually avoiding loops on the hexagonal lattice, and such loops are readily treated 
by CFT and SLE techniques. 

For general Q, the salient feature of Potts spin clusters is that the domain walls 
separating different clusters undergo branchings and crossings (see Fig. |TJ. These 
phenomena are however absent for Q = 2 with the above choice of lattice. It is precisely 
these branchings and crossings that make the application of exact techniques — such as 
CG mappings or Bethe ansatz diagonalisation — very difficult, if not impossible. The 
belief that spin clusters are indeed conformally invariant for other values of Q ^ 2 in 
the critical regime has even been challenged at times, but seems however well established 
by now [HI US]. 

Some progress has been accomplished in the Q = 3 case [UJ [TJ] by speculating 
that the spin clusters in the critical Potts model would be equivalent to FK clusters in 
the tricritical Potts model [151 E] ■ This equivalence has however not been proven, and 
is moreover restricted so far to the simplest geometrical questions [16] . The equivalence 
can also be understood as a relationship with the dilute O(n) model [T7] . 

The Potts model has been used recently to build a new class of 2D quantum lattice 
models that exhibit topological order [18J. Both FK clusters and domain walls between 
spin clusters are important in these models. We should also mention that spin domain 
walls have been studied numerically — and their difference from the FK clusters pointed 
out — in the context of the parafermionic models [T5] . 
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Apart from issues of branching and crossing, another major hurdle in the study of 
Potts spin clusters has come from the lack of a formulation that can be conveniently 
extended to Q a real variable. In the case of loops surrounding the FK clusters, 
this formulation led naturally to the introduction of powerful algebraic tools via the 
Temperley-Lieb (TL) algebra, and to the equivalence with the six-vertex model — the 
eventual key to the exact solution of the problem [5] . Factors of Q then appear naturally 
through a parameter in the TL algebra, or — via a geometrical construction — as complex 
vertex weights in the six-vertex formulation. Also for spin clusters can Q be promoted 
to an arbitrary variable: the weight of a set of spin clusters is simply the chromatic 
polynomial of the graph dual to the domain walls. From the point of view of the TL 
algebra, the domain walls are composite (spin-1) objects, hence more complicated than 
the FK clusters and loops. Recent work on the related Birman-Wenzl-Murakami (BWM) 
algebra [21] [20] suggests that this formulation might be amenable to the standard 
algebraic and Bethe ansatz techniques, although such a lofty goal has not been achieved 
so far. 

In this paper we report major progress towards the understanding of Potts spin 
clusters, both in the bulk and the boundary case. A brief account on the bulk case 
has recently appeared elsewhere [22]. Our results are of two kinds. On the one hand, 
we develop a transfer matrix technique which allows the formulation and numerical 
study of the spin clusters for all real Q. On the other hand, we identify the geometrical 
events that give rise to conformal correlations, and provide exact (albeit numerically 
determined) expressions for infinite families of critical exponents, similar to the familiar 
"L-legs" exponents j5j [7] for TL loops. Surprisingly, we find that geometrical properties 
of spin clusters encompass all integer indices (r, s) in the Kac table h T)S . An analytical 
derivation of our results appears for now beyond reach, in part because the algebraic 
properties of our transfer matrix are still ill understood. We do however provide some 
exact results based on an approach which does not involve a CG mapping, but rather 
the use of a massless scattering description. 

Apart from the bulk critical exponents [22], we report here the boundary critical 
exponents for free, fixed, and mixed boundary conditions. This means that the Potts 
spins on a segment of the boundary are restricted to take Qi values, where Qi can 
assume any real value and is in general different from the number of states Q taken 
by bulk spins. These results can be viewed as a further step in the programme [23] of 
classifying non-unitary boundary conditions in 2D geometrical models. 

2. Domain wall expansion 

The domain wall (DW) expansion of Q involves all possible configurations of domain 
walls that can be drawn on the dual of C (see Fig. [2|. A DW configuration is given by a 
graph G (not necessarily connected). The faces of G are the spin clusters. Since we do 
not specify the colour of each of these clusters, a DW configuration has to be weighted 
by the chromatic polynomial Xq(Q) °f the dual g ra P n G. Initially X&(Q) * s defined 
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Figure 2. A domain-wall configuration corresponding to a graph G, and its dual 
graph G. 

as the number of colourings of the vertices of the graph G, using colours {1,2,..., Q}, 
with the constraint that neighbouring vertices have different colours. This is indeed a 
polynomial in Q for any G, and so can be evaluated for any real Q (but X&(Q) * s integer 
only when Q is integer). For example, the chromatic polynomial of the graph G on 
Fig. § is Q(Q - l) 7 (Q - 2) 7 . The partition function Q can thus be written as a sum 
over all possible DW configurations 

rr NK ( -K\ length(G) , > 

G 

where N is the number of spins, and length(G) denotes the total length of the domain 
walls. 

The fundamental geometric object we consider is a connected part of a domain wall 
that separates two clusters. One can ask how the probability, that a certain number of 
such DW connect a small (in units of the lattice spacing) neighbourhood A to another 
small neighbourhood B, decays when the distance x between A and B increases. Each 
DW separates two spin clusters which connect A and B. There are in fact two types of 
such DW, depending on the relative colouring of the two clusters that are separated. If 
the two clusters have different colours, they can touch, so the DW is thin (see Fig. |3ja). 
If the two clusters have the same colour, then they cannot touch (otherwise they would 
not be distinct), so the DW has to be thick (see Fig. [3jb). 

Several different geometries are of interest. In the bulk case, the neighbourhoods 
A and B are at arbitrary, but widely separated, locations in the infinite plane. This 
is conformally equivalent to an infinitely long cylinder — a strip with periodic boundary 
conditions — with A and B situated at the two extremities. In the boundary case, the 
geometry is that of the upper half plane, with A being at the origin and B far away 
from the real axis. The boundary conditions are taken to be free, fixed or mixed on the 
positive real axis, and free on the negative real axis. In other words, the Potts spins 
can take Q\ (resp. Q) values on the positive (resp. negative) real axis. This geometry is 
illustrated in Fig. |4j it is conformally equivalent to an infinitely long strip, with A and 
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Figure 3. The two different types of DW, here shown in the bulk case (geometry of 
the infinite plane). A thin DW corresponds to the interface between two clusters of 
different colours (a), while for a thick DW the two clusters have the same colour. An 
illustration for the Q = 3 Potts model is given (left) as well as a schematic picture for 
non- integer Q (right). 

B situated at respectively the upper and lower extremities, and spins on the left (resp. 
right) rim taking Qi (resp. Q) values. We shall use the cylinder and strip geometries 
below to conduct our numerical calculations. 

We can now state the central claims of this paper. Consider the 2D Potts model 
for any real Q in the critical regime < Q < 4. Then the probability P that the two 
regions A and B, with separation i > 1, are connected by i\ thin DW and £2 thick 
DW decays algebraically. In the bulk case (geometry of the plane) the corresponding 
critical exponent is denoted h(Q, ^1,^2), and we have P oc x~ Ah ^' ei ' i2 \ Equivalently, 
on a long cylinder of size L x i with i ^> L, and A and B identified with the opposite 
ends of the cylinder, the decay is exponential: P oc q-^^I L ) KQA^) _ j n boundary 
case (geometry of the half plane) the critical exponent for free boundary conditions is 
denoted h(Q, £1, £ 2 )- 
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Figure 4. A thin DW in the half plane geometry for Q = 3. The left (resp. right) panel 
illustrates the case where the positive real axis supports free (resp. fixed) boundary 
conditions. 



Below we check these assertions numerically, and we observe that the numerical 
values of the exponents match the formulae 

h(Q,£i,£ 2 ) = he 1 ^e 2 ,2£ 1 , 

h(Q,£ 1 ,£ 2 ) = hx + 2(h-h),i+ih , (4) 
where we have used the Kac parametrisation of CFT 
(r-s k/A) 2 - (1 - k/A) 2 

K 

and 2 < k < 4 parametrises Q = 4 (cos G [0, 4]. 

It remains to give our results for the mixed boundary conditions with parameter 
Qi (of which fixed boundary conditions corresponds to the special case Q\ = 1). These 
follow by conformal fusion of the free-to-mixed boundary condition changing operator 
$ r0iSo that has been worked out in [23] with the operator ^i+2(£i-e 2 ),i+Ui that inserts 
the required number of DW with free boundary conditions. Details about this will be 
given below in section [5j 

We close this section with two remarks about the precise interpretation of the main 
result Q. 

First, in the expression Q for the bulk exponent h(Q, £±, £ 2 ), the case (£i,£ 2 ) = 
(0, 1) of a single DW (which must then necessarily be thick, due to the periodic boundary 
conditions) is special. In that case Q remains valid only if the DW is forbidden to 
wrap around the neighbourhoods A and B (in the geometry of the plane), or around 
the periodic direction (in the equivalent cylinder geometry). Without that restriction 
being imposed, we obtain another result: 

h(Q, 0, 1) = /io,i/2 (wrapping allowed) . (6) 

Second, in the boundary case the numbers £\ and £ 2 appearing in Q must be 
defined more carefully. Indeed, when there are £ = £\ + £ 2 propagating spin clusters, 
only the nature ('thin' or 'thick') or the £ — 1 DW separating them is clear by the 
above definition, whereas it is not yet clear how to characterise the 'half domain walls' 
that separate the two outermost spin clusters from the boundaries. The correct and 
unambiguous definition of £\ and £ 2 is as follows: Let the leftmost spin cluster contribute 



(5) 
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one unit to £\, and let each of the i — 1 subsequent spin clusters contribute one unit to 
l\ (resp. £2) if it has a different (resp. the same) colour as the cluster on its left. 



3. Transfer matrix formulation 

The DW expansion ^ may appear unwieldy and difficult to study numerically for non- 
integer Q. There nevertheless exists several Monte Carlo methods for studying the Potts 
model when Q is not an integer [2H 1213 EH], all of which are roughly speaking based on 
the FK cluster representation. It is possible to reconstruct the spin clusters from the 
Chayes-Machta algorithm [25] by performing the bond-adding step at zero temperature. 
Following our brief account [22] , this method has very recently been used [2Z] to verify 
a special case of Q corresponding to h(Q, 0, 1) = /i_i,o (with wrapping forbidden). 

As in [22] we shall however take a quite different route and resort to a transfer 
matrix construction. This has the advantage of linking up more easily with CFT [28J, 
and giving very precise numerical results for all the boundary conditions outlined above. 
Moreover, the ensuing transfer matrix formalism is no more complicated than the one 
[29] routinely used in the study of the FK clusters. 



3.1. State space 

Consider a strip of the square lattice of width L spins (boundary conditions will be 
detailed later). The basis states on which the row-to-row transfer matrix T acts contain 
one colour label c, per spin. By definition, one has Cj = Cj if and only if m = <Tj (i.e., 
the two spins on sites i and j have the same colour). The colour labels c, contain less 
information than the spin colours (7j themselves. For instance, any configuration in 
which the first and third spins have the same colour, no matter which one, and no other 
spins have identical colours, is represented by 

ci c 2 ci c 4 c L 

•••••••• (7) 

With these conventions we are always able to recognise whether to spins have the same 
colour or not, even if we do not know the precise value of this colour. But this is all 
that is needed to determine the Boltzmann weig hts in 0. 

For a row of L = 3 vertices, and for any (non-integer) Q, there are precisely five 
basis states 

C\C\C\ C\C 2 C\ C1C1C2 CiC 2 C 2 C1C2C3 



• • • . • • • . • • • . • • • . ••• 



(8) 



Note that the last state would carry zero weight for Q = 2, but apart from that the 
number of basis states for any given L will be finite and independent of Q. 

In general, the number of states is equal to the Bell number Bl of partitions of L 
objects, with exponential generating function 

00 „ 

£ TJ^ = ex P( eZ - !) • ( 9 ) 

L=0 
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Note that Bi differs from the Catalan number Cl of planar partitions, familiar from 
the FK cluster representation. In particular, {1,2, 1,2} is a valid state for L — 4. For 
L> 1 one has asymptotically Cl ~ 4 L , whereas Bl grows super-exponentially. When 
Q is integer, the number of states truncates to m Q L /Q\. 

Let us take the time evolution direction to be upwards, so that horizontal (resp. 
vertical) edges are space-like (resp. time-like). We can write T as a product of elementary 
transfer matrices, each represented symbolically as a rhombus surrounding a single 
lattice edge, and corresponds to the addition of that edge to the lattice. This edge 
links spins (shown as solid circles) on diametrically opposite sites of the rhombus. On 
an L = 4 square lattice with periodic boundary conditions this reads 




A rhombus corresponds to a vertical edge, and acts on a basis state s as follows. 
If exactly Q s distinct colour labels {c^} are used in s, then the new colour label d i of 
the spin Oi can be either unchanged = q (with weight e K ), or any one of the other 
labels already in use c[ = Ck (each with weight 1), or a new one c' { {c k } (with weight 
Q — Qs)- Note that this latter weight is in general non-integer, and is responsible for the 
correct computation of the chromatic polynomial X&{Q)- O n an example with Q s = 3 
this reads explicitly 



Cl c 2 c 3 



8 C „ c > e K • • • 

3 ' 3 CiC 2 C 3 



A. • • • + S C2 J • • • 

3 C 1 C 2 Ci 2 ' 3 CiC 2 C 2 



(l — £ci,c£ — ^c 2 ,c' 3 — ^3,4) (Q — Qs 



Cl C 2 C 3 



The last state could be written • • • , but this is equivalent to • » • in the basis 



A rhombus *C* adding a horizontal edge between vertices i and i + 1 corresponds 
simply to a diagonal matrix, with a weight e K if Cj = q+i, and 1 otherwise. 

With these rules at hand, one can write the periodic (cylinder geometry) L = 3 
transfer matrix for arbitrary Q in the basis 
hi ■ h 2 ■ h 3 ■ Vi ■ v 2 ■ v 3 with 

/ e K 



Q-l 



Vl 





e K 
1 





1 







\ 







Q-2 Q-2 



1 



+ Q 




as an instructive example: T 




1 
1 


3 A ' + Q 



3 7 
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and hi = diag(e , 1, e , 1, 1). The remaining matrices can be obtained from those given 
by cyclic permutations of the colour labels: 





/ 


e K 


1 













\ 






Q-l 


e K + Q- 


-2 












v 2 = 










e K 




1 


1 














1 




e K 


1 






V 








Q - 


2 


Q-2 


e K +Q 








e K 





l 










\ 









e K 







1 


1 




V3 = 




Q-l 


e 


K + Q- 


2 

















1 







e K 


1 






V 





Q-2 







Q-2 


e K + Q 


"3/ 


diag(e* l,l,e* 1), 


h 3 = diaj 


r(e K ,e K , 


1- 


1,1). 







(10) 



and h,2 

To get a non-periodic strip geometry instead of a periodic cylinder, one needs 
simply to omit the last horizontal edge which is responsible for the periodic boundary 
conditions. For instance, with L = 3 one would have T = h\ • h 2 ■ v i • v% ■ v%. This 
corresponds to free boundary conditions. Mixed boundary conditions are obtained by 
constraining the spins on the left boundary of the strip to belong to a subset of Q\ 
different colours, with < Q\ < Q. Spins on the right boundary remain free. For 
Qi integer this can be coded in the transfer matrix by considering that the number of 
colours in use in a given state is always at least Qi (so that in particular Q s > Q\). 
These first Qi colours are considered fixed, whereas other colours Q\ + 1, Qi + 2, . . . 
are defined only relative to the fixed ones, as described above. In particular, the choice 
Qi = 1 corresponds to fixed boundary conditions. 

The leading eigenvalue Ao of the transfer matrix T with periodic boundary 
conditions gives the ground state free energy f = — j- log A . This f coincides precisely 
with that of the usual FK transfer matrix p9j, even when Q is non-integer U Its finite- 
size corrections possess a universal L~ 2 term whose coefficient determines the central 
charge of the corresponding CFT 



3. 2. Correlation functions 

However, to obtain the desired two-point correlation functions of DW, we need to 
construct a variant transfer matrix T^^ 2 which imposes the propagation of l\ thin 
and £2 thick DW along the time direction of the cylinder (or strip). From its leading 
eigenvalue A^ 2 one can determine the energy gap A/^ 2 = — £ log(A^ 2 /A ) whose 
finite-size scaling in turn determines the critical exponents h(Q,£i,£ 2 ) and h(Q,£i,£ 2 ) 

[SHE]. 

To this end, the basis states ^ need to be endowed with some additional 
information about the connectivity of the spin clusters, ensuring the propagation of the 

I The reader may wish to check this fact on the L = 3 example given above. 
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desired number and types of DW. The crucial point is that we need to know whether 
two spins having the same colour label c$ also belong to the same cluster. Thus the 
states we use in the final transfer matrix have the form 



C\ C2 C\ C\ C5 C\ C2 C\ C\ C5 



meaning that two spins belong to the same spin cluster if and only if they are linked 
up in this pictorial representation. Of course, only spins with a common colour label 



can be linked up. Thus, in the left state of (11), the spins on vertices 3 and 4 are in 
the same cluster, but not in the same cluster as the spin 1. In the right state the spins 
1, 3 and 4 are all in the same cluster. These two states are different. In the transfer 
matrix evolution, each time two neighbour vertices correspond to the same colour, the 
corresponding clusters are linked up. 

Note that the possible ways of linking up the spins must respect planarity. For 
instance, in the L = 4 state with colour labels {1, 2, 1, 2}, one cannot simultaneously link 
the first and third spins, and the second and fourth, since this possibility is disallowed 
by planarity. 



To construct T^ 2 we modify the linked basis states (11) so that precisely i\ + £2 
spin clusters are marked. The colour labels of the marked clusters must respect the 
chosen values of £\ and £2- In order to conserve the marked clusters in the transfer 
matrix evolution, none of the marked clusters must be "left behind" , and two distinct 
marked clusters must not be allowed to link up. When a marked and an unmarked 
cluster link up, the result is a marked cluster. 

To summarise, the final transfer matrix thus keeps enough information, both about 
the mutual colouring of the sites and about the connectivity of the clusters, to give the 
correct Boltzmann weights to the different configurations, even for non-integer Q, and 
to follow the evolution of the boundary of a particular set of clusters. These boundaries 
are precisely the domain walls (see Fig. [3]) . 



4. Numerical results 



We have numerically diagonalised the transfer matrix in the DW representation for 
cylinders and strips of width up to L = 11 spins. We verified that the leading eigenvalue 
A in the ground state sector coincides with that of the FK transfer matrix, including 
for non-integer Q. As to the excitations A^ ^, we explored systematically all possible 
colouring combinations for up to 4 marked spin clusters, for a variety of values of the 
parameter k, and for several different boundary conditions (periodic, free, fixed, and 
mixed) . 

Finite-size approximations of the critical exponents h(L) and h(L) were extracted 
from the leading eigenvalue in each sector, using standard CFT results [281 [5], and 
fitting both for the universal corrections in L~ 2 and the non-universal L -4 term. These 
approximations were further extrapolated to the L — > 00 limit by fitting them to first 
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(2,0) 


(0,2) 


(3,0) 


(2,1) 


(0,3) 


(0,1)* 


2 


2.01(1) 


5.99(2) 


2.97(4) 




8.94(2) 


1.000(0) 


3 


4.01(1) 


7.99(2) 


6.02(2) 


8.04(2) 


11.98(3) 


1.500(1) 


4 


5.93(2) 


10.01(2) 


8.89(5) 


10.93(5) 


15.05(4) 


1.992(2) 


5 


7.77(4) 


12.09(8) 


11.6 (1) 


13.8 (1) 


18.2 (2) 


2.47 (2) 


Exact 


2(p-l) 


2(p + l) 


3(p-l) 


3p- 1 


3(p+l) 


p/2 



Table 1. Bulk critical exponents corresponding to five different DW configurations 
(•^1,^2), as functions of the parameter p — jz^r, along with the conjectured exact 
expression Q. The last column labelled (0, 1)* is for a single DW which is allowed to 
wrap around the periodic direction. The table entries give the value of |p|, when ^ 
is rewritten as h r _ s = {p 2 — l)/(4p(p + 1)), with error bars shown in parentheses. 





(2,0) 


(1,2) 


(2,1) 


(3,0) 


3 


7.00(3) 


19.0(2) 


15.16(4) 


11.13(4) 


4 


10.89(7) 


24.8(2) 


20.9 (1) 


16.74(5) 


5 


14.5 (1) 


31.1(2) 


26.7 (2) 


22.2 (3) 


Exact 


Ap - 5 


6p+ 1 


6p — 3 


6p — 7 



Table 2. Boundary critical exponents with free boundary conditions, corresponding 
to four different DW configurations (^1,^2), as functions of the parameter p = 
along with the conjectured exact expression Q. The table entries give the value of 
\p\, when ([5]) is rewritten as h r>s — (p 2 — l)/(4p(p + 1)), with error bars shown in 
parentheses. 



and second order polynomials in gradually excluding data points corresponding to 
the smallest L. Error bars were obtained by carefully comparing the consistency of the 
various extrapolations. 

Representative final results for the bulk case (periodic boundary conditions) are 
shown in Table [l] Recall that Q is only valid for (£±,£2) = (0,1) if the spin cluster 
is forbidden to wrap around the periodic direction. The sector where such wrapping is 
allowed is denoted (0, 1)* in Table [TJ and the corresponding critical exponent turns out 
to be h 0j i/2- 

Final results for the boundary case with free boundary conditions are shown in 
Table [2] The convention for the labels (£±,£2) is the one stated below Eq. (pi). 

5. Mixed boundary conditions 

We recall that mixed boundary conditions have been defined by imposing that the 
spins on the left boundary of the strip belong to a subset of Qi different colours, 
with < Qi < Q, while the spins on the right boundary remain free. In particular, 
Qi — 1 (resp. Qi = Q) describes fixed (resp. free) boundary conditions. As above, the 
convention for the labels (£±,£2) is the one stated below Eq. d6|. 
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An important point is that the fixed colour of the leftmost marked spin cluster may 
or may not be one of those allowed on the left boundary. For each choice of colours of 
the marked clusters, we define a sign e = 1 (resp. e = — 1) if the colour of the leftmost 
cluster is k\ < Q\ (resp. k\ > Qi), i.e., if the leftmost cluster is allowed to (resp. not 
allowed to) touch the left boundary. This sign closely parallels a construction in |23j, 
where the case e = 1 (resp. e = —1) was referred to as the blobbed (resp. unblobbed) 
sector. For free boundary conditions, only the choice e = 1 is meaningful (at least for 
< Q\ < Q, although it might be possible to give a meaningful definition of the model 
outside this range, by a suitable analytical continuation). The issue of the sign e = ±1 
is illustrated in Fig. |5j 




Figure 5. Mixed boundary conditions in the half plane geometry for Q = 3 and 
Qi = 2. The left (resp. right) panel illustrates the case e = 1 (resp. e = —1) where the 
propagating cluster, shown in blue colour, is allowed to (resp. not allowed to) touch 
the positive real axis. 



One would expect the mixed boundary conditions to be given by the conformal 
fusion of the operator $1+2(4 -^ 2 ),i+4^ a that inserts the required number of DW, and 
the free-to-mixed boundary condition changing operator $ ro ,s - The latter operator 
is responsible for the shift from Q to Qi allowed colours on the boundary. It has 
been worked out in [23] in the context of a different set of geometrical observables (TL 
loops), but since it should be representation-independent we can take it over here. The 
dominant exponents for mixed boundary conditions therefore follow from the standard 
CFT fusion rules [5] as 

(12) 

where h r0)So is the dimension of the free-to-mixed boundary condition changing operator 

The values of (r , s ) follow from [23], e.g., (—1, —2) for Qi = 1 and (— |, —1) for 
Qi = 2. We have verified these and other cases by explicit numerical calculations. In 



particular, the possible ambiguity of e versus — e in (12) is ruled out by the numerical 



results. Moreover (12) agrees with ^ for free boundary conditions (e = 1) as it should. 
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6. Massless scattering description 

While it is far from obvious to derive Q by CG methods, minor progress can be achieved 
using a rather different set of ideas. We restrict the discussion in this section to the 
bulk properties. 

We can learn about the dynamics of DW in the critical theory by using known 
information about the low-temperature (K > K c ) phase of the Potts model. Albeit 
non-integrable on the lattice, the corresponding deformation by the operator $21 is 
integrable in the continuum [30]. It can be described using a basic set of kinks K a b 
separating two vacua, i.e., ordered regions where the dominant value of the spin is a, 
resp. b. These kinks scatter with a known S-matrix related to the BWM algebra |21j. 
Importantly, the dynamics conserves the number of kinks: the process K a bKt> c — > K ac 
is forbidden (as in any elastic relativistic scattering theory), although kinks do appear 
as bound states in kink-kink processes. 

Many properties of these kinks can be calculated using integrability techniques. 
When the mass m — > (i.e., K — > K c ), the S-matrix provides a "massless scattering" 
description [3T] of some of the degrees of freedom of the critical theory itself. It is not 
entirely clear what a kink, which is well-defined for K 3> K c , becomes at K c , but it 
is natural to expect that thick DW are described by the propagation of two (or more) 
kinks such as K a bKb a - As for thin DW, the potential existence of regions where they 
are reduced to a single edge — that is a single kink — suggest they have to do with more 
complicated processes involving two kinks merging into one. We will not discuss them 
further here. 

It is an easy exercise to obtain the scaling dimension of thick DW using the massless 

scattering description. Indeed, the fact that the S-matrix satisfies relations from the 

(2) 

BWM algebra allows us to reexpress it in terms of the a 2 or Bullough-Dodd S-matrix 
[32], for which the thermodynamic Bethe ansatz was studied in |33j. This describes as 
well the dynamics of the field theory 

S = — [d 2 x \(d x §) 2 + (<9„$) 2 + g(2e'^ + e^*)! (13) 
8n J l J 

with f3 2 — j . Giving each kink the fugacity Q — 1 produces the correct central charge 

where each kink has a U(l) charge equal to 0, ±1. The scaling of the sector with charge 
j produces a gap Aj = 4-, so the leading dimension is Aj — (1 — c)/24 = hj/ 2 ,o- This 
agrees with /i_£ 2i o with £ 2 = j/2, so there are two kinks per thick DW. 
This argument validates Q for £1 = 0. 

7. Discussion 

The results Q and Q can be used to predict the fractal dimension of various geometrical 
objects related to spin clusters. It should be possible to observe these dimensions in 
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Monte Carlo simulations, and possible in real experiments. 



7.1. Fractal dimensions 

The dimension of a spin cluster follows from (|6| as 
d = min(2, 2 — 2/i 0)1 / 2 ) 

K e [2 ' i ] risi 

(8+k)(8+3k) r __ _ r8 xl l 10 ^ 1 



2 for K G [2, |] 

32K for«e[|,4] 



in agreement with [TUJ E] • 

The boundary of a spin cluster has dimension 



K 



4 = 2- 2/i_i >0 = 1 + - . (16) 

o 

This is an example of a duality relation: d^ for a spin cluster with SLE parameter 
k equals the dimension d PK of the boundary of an FK cluster at the dual parameter 
^* = ' ^ is moreover known that d^ K at parameter k coincides with the dimension 
d^ p of the FK cluster's external perimeter (with fjords filled in) at parameter k* [16J. 
Combining these, (16) means that d\> = d^ p , i.e., the dimension of the boundary of 



a spin cluster equals that of the external perimeter of an FK cluster, at the same k. 
This latter fact has very recently been verified numerically by Monte Carlo simulations, 
including for non-integer values of Q |27j. 

In the half-plane geometry, the intersection of the spin cluster containing the origin 
with the real axis has the dimension 

( 8 3/t\ 

d s = mm(l, 1 - h3,5) = mm l 1 ' 8- " - ^) > ( 17 ) 
as follows from Q. 



7.2. Thin and thick domain walls 

We have initially attached the epithets 'thin' and 'thick' to the two types of domain 
walls based on considerations in the microscopic model. Indeed, thin DW can narrow 
down to a single lattice edge on the dual lattice, separating spins of different colour. 
Since Q implies that the two types of DW scale differently in the continuum limit, one 
would expect that also the continuum geometrical objects can somehow be characterised 
as 'thin' and 'thick'. We now show that this is indeed the case. 

According to Q, the set of points where a thin DW has minimal width — i.e., one 
lattice spacing in the microscopic formulation, or the two clusters separated by the DW 
'come close' in the continuum limit in the sense of the small neighbourhoods defined 
previously — has dimension 

3 

dt = 2-2h 2A = — (4-k)(5«-4). (18) 
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This is analogous to the dimension of so-called 'red' or 'pivotal' bonds in the theory of 
percolation and of FK clusters. The corresponding result for a thick DW reads 

d 2 = max(0, 2 - 2/i_ 2 ,o) = . (19) 

The fact that d\ > means that thin D W are indeed thin, in the sense that they have a 
macroscopic number of loci of zero thickness in the continuum limit. Likewise, the fact 
that d 2 = means that the loci where thick DW have zero thickness is a set of measure 
zero. 

To summarise, there is therefore a consistency between the distinction of the two 
types of DW microscopically and in the continuum limit, and this ensures that they can 
(and do) scale differently. 

One can also obtain dimensions analogous to di and d 2 at the boundary. Indeed, 
define d\ (resp. d 2 ) to be the fractal dimension of the set which is the intersection of 
the real axis with the set of points where a thin (resp. thick) DW has minimal width. 
In other words, these are the dimensions of pivotal bonds on the boundary. By setting 
t\ = 2 and £ 2 = (resp. t\ = and £ 2 = 2) in Q we have then 

d± = min(l, max(0, 1 — ft.5,9)) 

for k e [2, f ] 
-8) for«G[f,3] (20) 
for k e [3,4] 




respectively 



d 2 = max(0, 1 - h- 3tl ) = . (21) 



7.3. The limit Q = A 

In the limit Q — > 4 (or k — » 4), we have d\ — > and d 2 — > 0. This means that 
the distinction between thin and thick DW disappears in that limit. This is of course 
consistent with the fact that k = 4 is a fixed point of the duality transformation k* = — , 
and so spin clusters and FK clusters have identical properties. 

A further manifestation of the indistinguishability of thin and thick DW is that the 
bulk exponent ^ becomes 

h - {il + 4)2 (22) 
Hx-l^Ux - t • \22) 

Indeed, this formula is invariant under the permutation of i\ and £ 2 . Moreover, the 
exponent for £ FK clusters reads C for Q = 4 [7] , and so a spin DW (whether thin 
or thick) behaves ClS db FK cluster, as expected. 

7.4- The case Q = 2 

For the Ising model Q = 2 (or k = 3) the absense of branchings means that one thick 
DW equals two thin DW. Indeed the bulk exponent becomes ([!]) becomes 

_ 4£ 2 - 1 

h ei ^i 2j2tl - — — — (23) 
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with i — i\ + 2^2 in that case, and this agrees with the exponent /i£/2,o for I loop strands 
[3 [7] in the dilute 0(1) model. 

7.5. The limit Q 1 

The limit Q — > 1 (or k — )■ |) is often studied in the context of FK clusters, since these 
become then bond percolation clusters. Setting Q = 1 for spin clusters has a more 
trivial meaning: all the spins are simply in the same state, Oi = 1. Therefore we should 



have d = 2 and d B for Q — > 1. This is indeed satisfied by (15) and (17). 

Apart from that, the Q — > 1 limit for spin clusters is more subtle. One can tackle it 
by considering configurations of bond percolation on the square lattice where one picks 
up one single percolation cluster. The spins on the vertices sitting on this cluster are 
given a certain colour. All the other spins in the system are given another single colour. 
Then the boundary of the coloured spin cluster is a single loop with fractal dimension 
df = 4/3, as expected from the duality k* = 16/ k. It is, however, remarkable that we 
can find such a simple construction of the spin cluster for Q = 1 on the lattice. For 
more interfaces one has to pick up more percolation clusters, and give them different 
colours. 

Finally, we note that for Q = 1 we have 

e 2 - i 

hh-ti&x = 24 , (24) 

which agrees with h 0t i/ 2 if one sets i = l\ + 3^2- This seems to indicate that for 
percolation, a thick DW equals three thin DWs (or TL loop strands) — a curious result 
for which we have no convincing explanation at present. 

8. Conclusion 

We have defined a set of geometrical observables based on the branching domain walls of 
the Potts model for any real < Q < 4. These observables are defined starting from the 
DW expansion ^ of the partition function. We have studied numerically these objects 
for non-integer Q using a transfer matrix formulation. Our results are compatible with 
conformal invariance of these observables and we have given a set of conjectures for the 
corresponding bulk and boundary critical exponents. 
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